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SUMMARY 

This report develops results pertinent to predicting the 
performance of convolutionally encoded binary phase- shift keyed 
communication links. The details of the development are 
provided in foux' sections. These sections are concerned with 
developing the bit error probability performance degradations 
due to PN despreading by a time -shared delay locked loop, the 
Costas demodulation process, symbol synchronization effects and 
cycle slipping phenomena in the Costas^oop. In addition, Costas 
cycle slipping probabilities are studied as functions of Doppler 
count time and signal-to-noise conditions. The effect of cycle 
slipping in the symbol synchronizer is also studied as a function 
of channel Doppler and other frequency uncertainties. 

Section one develops the model of a suppressed carrier tracking, 
spread spectx'um receiver/transponder in which the received 
signal Is a convolutional/Manchester encoded binary data stream 
that is spi'ead by a PN code to produce a binary phase-shift keyed 
(BPSK) carrier. The combined signal-to-noise degradations of 
the PN despreader^s time -shared delay locked loop and the Costas 
phase-locked loop are developed as fimetions of the various 
system parameters and channel conditions. In addition, the theory 
of cycle slipping in the Costas loop is developed in terms of the 
channel Doppler dynamics, the additive noise, the minimum 
squaring loss, the PN code tracking loop jitter and the duration 
of the Doppler count normalized by the link data rate. Cycle 
slipping probabilities are developed as functions of the energy- 
per bit to noise ratio, the nmnber of cycle slips, the Doppler 
count time normalized by the link data rate and the receiver's 
static phase error. Various cycle slipping probabilities ar e 
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grapHcally demonstrated as a function of the Doppler coxint time. 

Section two develops the bit error probability expression that 
results from the noisy reference loss due to jitter in the PN 
despreader's time-shared delay-locked loop and the noisy refer- 
ence generated by the Costas demodulator. For the K=;7 rate 
one-third convolutional code, bit error probability curves are 
graphically displayed as functions of the data rate, static phase 
error and energy-per-bit to noise ratio. Dink performance 
degradations ax*e easily assessed from the curves; in addition, 
any luiown losses due to filtering of the signals in the receiver 
can be easily added to the loss indicated by the curves presented. 

Section three of this report is concerned with symbol synchro- 
nization in the Network Signal Processor when Manchester coding 
is used. New and fundamental results pertaining to symbol 
synchronization for Manchester coded messages are derived. 

The symbol synchronizer structvires are found to he functions 
of the signal-to -noise ratio and data transition density. The 
minimum rms symbol sync jitter is determined, computed and 
comparisons are made for various phase detector characteristics. 
A general theory for cycle slipping in symbol synchronizers 
is given; it is shown that under most cases of interest, the 
effect of slipping in the symbol synchronizer does not give 
rise to a significant performance degradation. This section 
is also used to study the effects of synchronizer stress (due to 
clock uncertainties and channel Doppler) on bit error probability 
performance. 

Section four of this report evaluates the degradations on bit 
error probability perforniance due to the noisy clock reference 
created by the symbol synchronizer. Particular graphical 
results are presented for the special case of a K=;7 rate one- 
third convolutional code. The effects of data transition density 



on symbol synclii‘onizer/decoder performance is also studied and 
graphical results for cases of greatest interest are presented. 




1.0 SUPPRESSED CARRIER TRACKING FOR SPREAD SPECTRUM 

RECEIVERS OF MANCHESTER CODED MODULATION 

1.1 Inti'oduction 

Various communication systems transmit information in the 
form s(t) = m{t)Sp^(t)siii[i«Qt t 0(t)^ where Spj^(t) is a PN 
modulation and m{fc) is a linear modulation (coded or uncoded) 
which possesses no dc component in its power spectrum. Since 
no residual spectral component esdsts at the carrier frequency 
it is not possible to use the conventional phase-locked loop for 
establishing the required coherent reference. 

A number of methods have been proposed for generating a 
reference car’rier from the received waveform which contains a 
suppressed carrier signal component. Of these, the three most 
popular are the squaring loop, method (Refs, 1-5), the Costas loop 
method (Refs, 1,Z, 4-6), which is theoretically equivalent to the 
squaring loop, and the decision-directed feedback loop (Refs. 1,2, 
7-9) which first estimates the modxjlation, and then uses this 
estimate to eliminate the modulation from the carrier, leaving, as 
nearly as possible, an momodulated sinusoid which can be tracked 
with a phase-locked loop. Each of these methods has advantages 
and disadvantages (practical or otherwise) with respect to the 
other two. What ultimately determines which method is used in 
a particular application is a tradeoff between the requirements 
on performance capability and the ease and cost of impiementation. 

As is well-loiown from previously published documents (Refs, 
10-12) suppressed carrier modulation is employed on both the 
TDRSS-to-Orbiter (Forward) and Orbiter-to-TDRSS (Return) 

S -band links. In the case of the forward link, the convolutionally 
encoded, time-division multipleiced (TDM) data sequence is 
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convex'ted from KRZ-Li to Bi-phase-L foi’mat, modulo-2 added to 
a PN NEZ--Lt code sequence, and the resulting spread-spectrum 
sequence used to biphase modulate the transmitted carrier* After 
transmission through the TDRSS, the Orbiter receiver first 
despreads the signal by acquiring and tracking the PN code. 
Following this procedure, the baseband data symbol sequence 
is restored by demodulating the suppressed carrier signal with a 
Costas loop. In the return link, similar carrier modulation and 
demodulation processes taJce place with the exception that the 
signal is neither PN spread nor despread. 

It is clear from the above that in order to assess the degrad- 
ing effects of the carrier regeneration process on the overall 
system performance (including the average error probability 
performance of the data decoder), a thorough understanding of 
the acquisition and tracking performance of Costas loops and 
their optimum design for use in the Shuttle relay linlc environment 
is x’equired, 

^ Performance Analysis of the Costas Lioop in the Suppressed 

Carrier Mode - Es^ect of Spread Specfcrum Modulation 

‘1.2,1 Motivation 

The tracking performance of a Costas loop in response to a 
biphase modulated suppressed carrier input is well-documented 
in the literature (Refs, 1,2, 4-6), All of these analysis have 
assumed that the in-phase and quadrature arm filters have 
sufficiently wide bandwidths so as to pass the data modulation 
undistorted. In practice, the bandwidths of these filters are 
more commonly chosen on the order of the data rate and thus the 
above assumption is strictly speaking invalid. In Refs, 13 and 14, 
the effect of arm filter distortion has been studied in detail. 
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graphically demonstrated, and compared with the results given in 
the literature which have neglected this important effect. In fact, 
careftil control of the distortion term in any design gives rise to 
the highest noise immunity achievable with passive arm filters. 
Even these analyses, however, have not included the degrading 
effect of a spread spectrum modulation on the signal being tracked. 
It is the purpose of this section to augment the previous work 
on Costas loops hy including the effect of a PN spread spectrum 
modulation on tracking performance. The approach will parallel 
the development taken in Ref, 13 and thus much of the detail 
given there will he omitied in the presentation here. Finally, the 
results obtained will be adapted to cover the case where the 
Manchester coded data is also convolutionaUy encoded, 

1.2,2 C-ostas Lioop Model and Analysis 

Consider the transponder illustrated in Fig, 1-1 whex'e the 
Costas loop portion is enclosed v/ithin the dashed outline. Since 
our main interest here is in the performance of the Costae loop 
itself, it is sufficient for us to model the signal at its input and 
concentrate on how the loop processes this signal. Referring to 
Fig. 1-1, if the received signal Xj^{t) at point @ is modeled as 

“ /2s's^^(t)m(t) sin ®{t) + n.(t) (1) 

then, using straightforward trigonometric manipulations, one 
finds that the signal at the input to the Costas loop (point 0 ) is 
given hy 

5(t)] + 

( 2 ) 

In Eq. (1), 5(t) suJq t -[■ 0(t) with OU the radian carrier fi’equency 




















and 0{fe) ^ + 6 q the input phase to be estimated, m{t) is the data 

modulation (a ±i digital waveform), s (t) is the received PN 

jprv 

modulation, and m(t) is the additive channel noise which can be 
expressed in the form of a narrowband process about the actual 
frequency of the input observed data, i. e, , 


n.{t) = /T{N (t) cos §{t) - ISf (t) sin 9(t)] 
X c s 


(3) 


where N (t) and N (t) are approximately statistically independent, 

C •• ■ s 

stationary, white Gaussian noise processes with single-sided 
noise spectral density w/Ha (see Ref. 2) and two-sided band- 
width B < WJ /^TT. In Eq, (2), s_ (t+T ) is the PN reference signal 

U U JT S ^ 

derived from the PN tracking loop, is an equi- 

valent loop gain, and 

n^(fc) “ v^"{N^(t)cos[f(t) l(t)] - Ng{t)sin[*(t) ^(t)]3 (4) 

Also assumed in Eq. (2) is that the bandpass filter which precedes 
the Costas loop is sufficiently wide as to pass the data modulation 
m(t) undistorted. Denoting the in-phase and quadrature detector 
(multiplier) gains by then the output e^(f^) tiie quadrature 

phase detector is (ignoring second harmonic terms) 

A 

cos H) = K^KJKJ^/S’Rp^^(T^)m(t)-N^(t)] sin cp(t) 

+ KK,K N' (t) cos cp(t) (5) 

e 1 m c 

while the in-phase arm phase detector output is 

A 

^s'^> = sin = K^KjK^[7S'K.pj^(T^)m(t)-N^(t)]cos,¥(t) 

-K K,K N*(t) sincp(c) (6) 

e i m c . 





where q3(t) =$(fc) - 
A 


(M-f 1 )N+1 ^ 

-- ■ S(t) is the loop phase e«or. 


MN 


= Sp^{t)Spj^(t+T^) is tlie autocorrelation ftuictlon of tiie PN process 
which, is approximately given by 




1 . 

A ' 


0 


iT^l^A 


elsewhere 


(7) 


where A denotes the PN chip interval, and N' {t^ and N* (t) are 

c s 

equivalent noise processes defined by 

n;( 9 = Spjj(t+T^)N^(t) 




(8) 


In writing Eqs* (5) and (6), we have again assumed that the data 

modulation is passed undistorted by the bandpass filter following 

point @ . Furthermore, we have assumed {and less justifiably) 

that this filter has no effect on the PN autocorrelation function. 

In practice, the autocorrelation function of (7) is somewhat spread 

by this filter and furthermore, the noise processes in (8) should 

reflect the filtering effect on the PN demodulation signal s„..(t+T ). 

PN e 

Without producing undue complexity into the analysis, we shall, 
to a first order approximation, ignore these effects. 

After low -pass filtering with the in-phase and quadrature - 
phase arm filters, the phase detector output signals become, 
respectively, 

^ N'‘(t) 

m(t) s 

" ^e^l^m^^ ^PN^'^e^ G{p)m.(t) - G(p)N^(t)] sin cp(t) 

N"(t) 

~ (9) 
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// 


m(fc) 


N“(t) 

s 


• Z.Jt) = K^KjK^Cyg'R^„(T^)G(p)m(t) - G(p)N;^(t)] cos (p(t) 


•PN' e 


S 

N” (t) 
c 


- K K, K [G(p)N‘ (fc)l sin cp (t) 
e 1 m c 


( 10 ) 


Here the signal in(t) denotes the signal emerging from passage of 
the data modnlation m(t) through the low -pass arm filter G(p)j, and 
likewise N^(t) and N”{t) are the equivalently filtered versions of 
N* (t) and n’ (t)^ r espectively* Multiplying the two low -pass filter 

C 5 

outpvits (assuming this multiplier has unit gain) gives the dynamic 
error signal 

2 2 2 

K K K , 

(t)}sin 2(p(t) 


+ K^K^K^N^(t)L'S Rpj^(T^)ii(t)-N^'{t)3 cos 2(p(t) 


( 11 ) 


The instantaneous frequency of the VCO output is related to 2 ^(t) by 

^ = K^[F(p)SQ(t)] + Uj, (12) 

and hence the stochastic integro -differential equation of operation 
of the Costas loop becomes 


2 =: 2Q^ - KF(p){SRp^,(T^)m^(t) sin 2cp(t)+ v^ft, 2tp{t)]} (13) 


dt 


2 ''2 


A 2 2 2 

where K = K K, K K„, and 
e 1 m V 


v^[t, 2«p(t}]^[-N”^(t)+Nf (t) ~2VSll^^(T^)m(t)N|;(t)]sin2cp(t) 


+ [2VSRp^(T^)m(t)N”(t) -2N"(t)N“(t)]cos2,:p(t) 


(14) 


In arriving at (14), we have made the practical assumption that 


- 10 . 


/?■ .f? 
'GiLiiriL^otn 


/) f? 


the data rate =1/T is large relative to the equivalent loop 

s iS 2 ^ 

bandvnidth W , and thus m (t) can be replaced by its mean- 


squared value, io e. , 


■^(t) 6 J S^(f)|G(jETF£)fdf 


with S^{f) denoting the power spectral density of the data modula- 
tion m(t) ... 


, 1.3 Steady-State Tracking Performance 


. Usmg Fokker -Planck techiiiques, the steady -state probability 
density function (p.d. f. ) p(2^) of the modulo Z-n reduced phase 
error Z^ can be determined from (13). Assuming a loop filter of 


the form 


F(s) = 


I+ST2 


F, = - 

1 r- 


then, (Ref, 1, Chap. 2 and Ref. 2, Chaps, 9, 10) 

2^+2tt 

p(2^) = Cq exp02{j) -t a cos 2^)J e3cp(-p2x-a cos 2x:)ds: (l7) 

■ . 24 ... 


where Gq is a normalisation constant and 
“ = > - ~2- 

I 

P = ■ ®I^pi^(Tg)ni^{t)K(l-Fj^)sin^]+asijrz^ 

.hi 


with 


2 A 2 ■ 

r = second order loop damping parameter = SR^.It )m (t)KF, t. 

o Jt^JN 0. 1 j 

= 4C^ 





f"' 

)] 


. p. 


C = loop damping 

p> = effective signal-to-noise ratio in the loop bandwidth = 

p =: equivalent signal-to-noise ratio in the loop bandwidth of 

second -o rder PLL = 2SR^ _(t )/N_'W 
^ 3 U Xj 

a = sin 2<j) ^(sin2^) 


=: loop squaring loss = 


N 




N 


sq 


0 ) 

= zj (T)dT 


(19) 


oa Z 


and 


(t) = V2(t, 2cj))v2(t+T, 2(|) = 4[SRp^(T^)R^(T)Rj^ri(T)+R^ri(T)] (20) 

2 c c 


where 


N„ « 


and 


Rj^,.('r)- N^(t)N”(t+T) =~ J |G(j2rrf)lV^^% (21} 

c -c» 

R^(t) = ^{t)m(t+T) =: J S^(f)|G(j2TTf)fe^^^^'^df (22) 

■-13 


In arriving at (20) we have made use of the previous assumption 
that the arm filters are narrowband relative to the input bandpass 
(IF) filter. The probability density function in (l7) will be useful 
in what follows. 

The squaring loss el can be derived in terms of basic system 

Xj 

parameters^ Substituting (21) and (22) into (20) and using the 
definition of in ( 19 ) giw 
the squaring loss formula 


definition of in ( 19 ) gives^ after considerable manipulation, 
Xj 


eepboducibility of the 
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1 = 


B./-? 

V . 1 ^ 

^ 2 R‘D 

s m. 


^ A 2 

where we have used D =m (fc) to denote the modulation distortion factor, 
2 

R^'=SRp^(t^)T^/Nq =R^Rp^{T^) is the effective data (symbol) signal - 

to-noise ratio, B. denotes the two-sided noise bandwidth of the 
1 

arm filter G(j2TTf), i, e.. 


B k r |G(j2Tlf)|^<if 


K Is a constant dependent only on the filter type and is defined by 
Xj 


L £3 


09 A 

J 1 g( 32 ti£)| df 

J lG(j- if 


and is a constant dependent on both the baseband data power 
spectrum and the filter type, i.e. 


D ® 


^ • A 

I S^(f)lG(}2nf)| df 

-CO 

CO 

J S^(f) |G{j 2 Tlf)l df 


Typical values of K_ for well-known filter types may be found in 
• . Xj 

Table 2-1 of Ref, 1, Chap, 2, For example, K = 1 for an ideal 

1-f 

low -pass filter while K s (2n-l)/2n for an n -pole Butterworth 
filter. Since, via Eqs, (15) and (26), the modulation distortion 
factor and the constant Kj^ respectively depend on the 
baseband data power spectrum S^(f), the format of the baseband 
data encoding must be Specified before these quantities can be 




comp-uted. The case of inherest here is when the modtilation m(t) 
is a Manchester coding of equiprobable, independent transmitted 
symbols. The power spectral density S^{f) for such a data 
modulation is, Ref. 1, Chap* Z, 

S^(f) stn^nfT /2) 

rp " 2 

. (TTfT /2) 
s 

Recalling that an n-pole Butterworth filter is characterized by 
the transfer function 


lG(j2trf)l = 


1+Cf/f ) 


2h 


(i 


where f , the 3 dB bandwidth, is related to the two-sided noise 
c 


bandwidth B. of the filter by 
1 



then, substituting (Zt) and (28) into (15), (25), and (26), the 

modulation distortion factor and constants K„ . and K_ can be 

m L» D 

computed by numerical integration as functions of the ratio 
B,//? . Using these results, Fig, 1-2 illustrates the squaring 

X 5 

loss ^ vs B,/P with R* as a parameter for a two-pole 
L I s s 

Butterworth filter, i,e* , n = 2, We observe that for a fixed 
R^ there exists an optimiim noise bandwidth B, for the arm 
filters in the sense of minimizing the squaring loss* These 
values of optimum arm filter bandwidth occur in the vicinity of 
the Nyquist bandwidth. 

In designing carrier tracking loox^s, the most commonly 
used measure of performance is mean-sqiiared tracking jitter. 
This quantity is ordinarily defined as the variance of the 
probability density function of the loopfs modulo Ztt reduced error 


/I 
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and for large loop signal-to ’noise ratios becomes equal to the 
reciprocal of the effective loop signal- to -noise ratio. In the case 
of interest here, the tracking performance of the carrier -tracking 
loop depends, in addition, upon the timing error in the PN track- 
ing loop which itself has a statistical characterization. Thus, 
the mean-squared tracking jitter for the Costas loop of Fig, 1-1 
is given by 


TT 1 , 

0^1 ~ J J (2{(i)\{24)p (T^/A}d(2<|.)d{T^/A) 


(30) 


where p(2^>) is given by (17) and p (t /A) is the probability density 
function of the normalized timing error T^/A in the PN tracking 
loop. Although p^(t /a), in general, depends upon the specific 
PN tracking loop imxJlementation, e. g,, a noncoherent delay- 
locked loop or T 'dither loop, for large PN loop signal -to-noise 
ratios, p^(r^/A ) is approximately Gaussian distributed, i.e. 


p (t./A) “ 
T e 



exp 


(t /a) 
© 


2 -1 


2o 


(31) 


and is thus characterized by a single parameter the standard 

deviation of r /A. Furthermore, for large Costas loop signal-to- 

■. ' © 

noise ratios, 


^2c|. 


_L 


(32) 


where 


P' “ f P'p(T /A)d(T /A) 
J ^ T e e 


(33) 


Substituting (7), (19), (23) and (31) into (33) results in an expression 
for the average effective signal-to-noise ratio in the Costas loop 
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[s 0 ^ /■ 







D 


m 


Kfi+Kl, 


B./^‘ 

1 s 




7 


exp 


'zna 


2 n 


n 


Ea 


T -J 


dT 


n 


(34) 


where we have introduced the notation t = t /A and recognized the 

n e 

fact that the symbol signai-to-noise ratio R can be written as 

5 


R - 
a 


2S 


N„W_ 

^ 0 L 


(35) 


A ^^s 

with 6 = ", '' ' the ratio of the symbol rate to single-sided Gostas 

s W 

loop bandwidtii. The value of B.//? to be used in (34) is that which 

^ IS 

minimizes (as in Fig, 1-2) for a given R . Since r’ = R 
Xj s s s 

when T =0, these optimum values oi^jR can be determined 
e is 

from Fig. 1-2 and remain constant over the integration on in 
(34). The values of D and to be used in (34) are also 
evaluated at these optimum values of Figs, 1-3 and 1-4 

plotp* vs fitted and a parameter. The dotted 

curve in each of these figures corresponds to the case of no 
squaring loss or PN degradation, i.e., p* = (2 S/NqW^)/ 4. Thus, 
the distance between this dotted curve and the curve labeled 

a = 0 represents the si gnal-to -noise i*atlo degradation due solely 

. ■ T 

to the Gostas loop squaring loss. The remaining curves represent 
the additional degradation due to imperfect PN tracking. These 
curves are very useful in assessing the combined performance of 




the PN despreader and the Costas loop. 

1,4 Characterization of the Cycle “Slipping Probabilities 

In this section we shall characterize the probahility that the 
2^ process in the Costas loop slips one or more cycles in an 
Interval T seconds. To do this, we shall as snme that the cycle 
slips which tend to increase and decrease 2cp by 2tt form independ- 
ent Poisson processes with rates of oecui'rence N , and N 
respectively. Then the probability tliat a net mxmher, say N = n, 
of cycle slips occur in T seconds is given by (Ref, 2, Chap* 9) 

P(N=n) = exp{-ST)I (2Tv^Tir) (36) 


N 


for n = 0, ±1, ±2, . , , where S =j]Sf , + N is the total number of cycles 
slipped per unit time independent of direction and I^(x) is the 
modified Bessel function of order n and argviment The result 
(36) is shown by forming the convolution of the difference of two 
independent Poisson processes. The fact that the p. d.f, for 
cycle -slipping can be approximated by a Poisson process has been 
observed experimentally. 

To rewrite Eq, (36) in terms of the system parameters 
previously defined we note from Ref, 2 that 


!!h 

N 


~ exp(2TT0) 


W. 




exp(j:Trg) 

2nV J 


_ W 

C? X 4. 

s = -p r 


coshnp 

TT^Uj3(p‘)f 


(37) 


^'We shall for simplicity omit the effect of PN degradation in these 
calculations. Also a first-order loop will be assumed. 




Thus, from (37) we get 


W 


TipT" = — — ^ ^ 

+ - E/p'|yp')r 


(38) 


Substituting (37) and (38) into (36) and making use of the definition 

‘bf 6 following (35) gives the desired result 
s 


P(N=ni = exp 


/ T \ 

f 

\ s/ 




= exp 


UTTp - 


Z cosh TT^ 

P'TV^6jl.p(P')f 


/ t\ 

It / 

\ s / 


T 


(39) 




I ,2 

Various methods for evaluating |I (p')l numerically are given in 

JP 

Appendices n and IV of Ghap^ 9 in Ref. 2, For example, 

= SejipVpn) J 


For a first-order loop, operating with steady-state phase error 

<1» , we have that 

^ss 


P = p' sin 


ns 


(41) 


Also, from (19) the effective loop signal-to -noise ratio p* can be 

expressed in terms of the symbol signal- to-noise ratio R = 

s 

ST /N. by 
s 0 ■' 


R 6 

p? _ — ®_® j> 

H - 4 


( 42 ) 


or for convolutionally- encoded data of rate 1 /v 




p« = 


K6 

b s 

4v 


(43) 


Thus/ substituting (4X) and (43) into (39), one obtains an ejipression 

for the probability of N aycle slips in T seconds as a function of 6 , 

s 

T/T , R (or R, ), and ^ • Figs, 1-.5 to l-lO plot P(N"n) vs T/T 

for 6 and i|) fixed and R as a parameter. Values of n ~ 1 and 
n - 2 are chosen as representative of the physical situation since 
the probability of three or more cycle slips is insignificant relative 
to the other two probabilities. Also, the rate 1/3 Shuttle code is 
assumed and R^ is allowed to range over the values -1 dB to 4 dB, 
The value used in aU cases is the optimum as determined 

from Fig. 1-2, 
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2. 0 AVERAGE ERROR PROBABILITY PERFORMANCE OF 

CONVOLUTIONALLY CODED SYSTEMS IN THE PRESENCE 

OF IMPERFECT CARRIER SYNC 

Z, 1 lL?.troductiQn 

Having thns characterised two sources of imperfect synchro- 
nization that degrade system performance what remaina is to 
evaluate their effect on the average error probability performance 
of the data decoder which in the case of the Shuttle forward and 
return S-hand relay links is a maximum-likelihood (Viterbi) 
decoder of a rate one -third, constraint length seven convolutional 
code. A general solution of this problem is indeed quite difficult; 
however, depending on the relation among particular system 
operating parameters, e.g, , symbol rate, decoder memory, loop 
bandwidths, etc,, there exist approaches which yield veritable 
results in the face of the underlying assumptions and approximations 
made. 

One approach which has been taken by several avithors {Refs. 
1,2) with regard to evaluating performance degradation due to 
imperfect carrier synchronization is as follows: (a) determine 
(by simulation studies) the ideal bit error probability performance 
curve for the particular convolutional code under consideration, 

(b) curve fit (using perhaps a least squares approach) a mathe- 
matical function to the simulation data, (c) degrade the signal-to- 
noise ratio in this function by the cosine squared of the carrier 
synchroniza.tion phase error, and (d) average the resulting 
function over the p. d.f, ,of the phase error as determined from 
the analyses performed on the carrier synchrcnization loop. 

Inherent in this approach, is the assumption that the loop phase 
error is approximately constant over the number of symbols 
required to arrive at a single bit decision, i.e, , the Viterbi 
decoder memory length in symbols. Typically, this memory 




runs on the order of five or ste constraint lengths for which a 
constraint length seven code is about forty symbols. Thus a 
symbol rate of 216 K^sps would require that the two-sided carrier 
tracking loop bandwidth, he much less {e.g., a factor of 

100 or more) than 216/40 = 5.4 KHz. 

We begin by finding a suitable curve fit model for the 
Shuttle K = 7, rate 1/3 convolutional code. This code, which is 
the ”best’' rate 1/3, constraint length 7 convolutional code (Ref, 3), 
is cliaracteriaed by the generator matrix. 


G “ 


L 


11110 0 1 
1 1 0 0 1 0 1 
1 0 1 1 0 1 1 


( 1 ) 


and is the identical code used by the LINKABIT li‘V70l5LR 

Viterbi decoder which employs 3 -bit soft quantisation of the 

received data. Using previously published simulation results 

(Ref. 4) we find that over the range of bit signal-to-noise ratios 

from -1 dB to 4 dB (corresponding to a range of bit error 

-i -5 

probabilities from approximately 2x10 to 10 ), a- second- 

order polynomial fit is quite accurate, i.e. 


where the constants a^, and a^ are chosen to satisfy three 
data points on the simulated curve. In particular, simultaneous 
solution of the three resulting linear equations gives 


= -.1206 


^2 “ "* 


5 


(3) 


Using this curve fit model, numerical results can be obtained for 
the noisy x-eference effect of the Costas loop on the bit ei*ror 
probability performance of this particular code. The results will 
be plotted with the ratio of bit rate to loop bandwidth held constant 
and static phase error (caused by Doppler offset) as a parameter. 

In all cases, the Costas loop is optixuized in the sense of minimizing 
the loop's squai'ing loss (or equivalently, its mean-squared jitter) ' 
by proper choice for the arm filter bandwidth relative to the symbol 
rate. A two-pple Butterworth is chosen for the arm filters, 
representing a reasonable compi*omise between filter complexity 
and modulation distortion. 

2,2 An Analytical Evaluation of Performance Degradation Due to a 
Noisy Carrier Sync Reference 

Based upon the comments made in the introduction, we can 
immediately evaluate the average bit error probability performance 
in the presence of a noisy carrier sync i*e£crence from the following 
numerical integration: 

TL 

2 2 

TT 

"2 


where from (2) and (3) 



5(l01og^QR|^ 


2 

cos 


-1.0706 


(S) 


and p(2^ ) is the p. d.f. of the suppressed carrier -tracking loop 
phase error. Implicit in (4) is the fact that the 180*^ phase 
ambiguity, which is cliaracteristic of all suppressed-carrlei* 
tracking loops, has been perfectly resolved. Techniques for 
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accomplishing this resolution ai'e discussed in the section on 
average symbol transition density of convolutional codes. 

When a Costas loop is employed for establishing the required 
carrier sync reference, then the p. d.f. p(2{|> ) is given by 

p(2(^) = exp^‘ cos f s3J:p[-p* cos Ey-(3{Ey)3dy; 

(6) 


where Cj^ is a normalization constant such that 


and 


Tf 

J p( 24 )d^ ” p(^4)d4 “ 1 

-TT 


, 2S 

P " — = NlW 


0 L 


J3 = p ’ sin (|) 


ss 


(7) 


( 8 ) 


In (8), the parameter is the loop’s squaring loss, W is its 

Li Xj 

two-sided noise bandwidth, and c[f is the static phase error 

S 3 

caused by Doppler offset. Defining the ratio of symbol rate to 
single-sided loop bandwidth by 


2 /?, 


6 r. , 

s w 


:b 


(9) 


then 


P = 


R 5 

b 5 
V 


(10) 


where 1/v is the rate of the code. 

Figures 2-1 and 2-2 represent plots of bit error* probability 
vs energy per bit to noise ratio. These results were 


iOBUCIBEirrY OF THE 
IS POOR 


- 34 - 







/I 

obtained via numei'lcal integration of (4) nsing (5) and (6), These 
curves are illustrated for the Shuttle K = 7 and v = 3 convolutional 
code* For 6 = 432 or 216, notice the degrading effects due to the 

B 

static phase erroi'. In addition, it should be noted that other 
degradations due to filters in the signal path can easily be 
incorporated and Figures 2-1 and 2-2 used to indicate the total 
loss due to the noisy carrier reference, when the various filters 
in the receiver are assumed to be presents 
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3, 0 SYMBOL SYKCHRONIZATION OF MAFfCHESTER CODED 

DATA 

3. 1 Introduction 

The technical literature ahoimds Vi^ith information relating to 
the problem of establishing symbol synchronization in binary 
baseband digital communication systems transmitting NRZ-L 
data (Refs, 1-15), By comparison, very little has been said 
(Refs* 16-18) regarding the extraction of symbol sync from Bi- 
phase-L (Manchester coded) data. In many applications (e. g. 
Shuttle relay links), however, the advantage of an increased 
average transition density (a factor of three for random data) 
offered by Manchester coding overslmdows its potential dis- 
advantages such as the introduction of a half-symbol ambiguity 
requiring resolution and the approximately doubled bandwidth 
occupancy. Thus, in these instances, the problem of providing 
symbol synchronization in the receivers of such systems requires 
special attention. 

3. Z Maximum A Posteriori (MAP) Estimation of Symbol Sync 

Regardless of the data format, the question of what is the 
"optimum” symbol synchronizer is a natural one to ask* The 
answer to this question, of course, depends highly upon the 
criterion of optimality one selects. The maximum a posteriori 
(MAP) estimation criterion represents one approach to the problem. 
Although this approach Tinfortunately avoids a direct attack on the 
important problems associated with symbol sync acquisition time, 
acquisition range, and maximum probability of acquisition, it does, 
however, afford a great deal of insight into arriving at closed 
loop synchronizer configurations whose high and low-signal-to- 
noise ratio approximations are physically and practically 





om. 




3fealizaT3le, and yield good, if not the best, tracking performance. 
In fact, we shall demonstrate that, under certain practical con- 
ditions^ the rules governing the design of the symbol synchronizer 
motivated by the MAP estimation approach are indeed identical 
to those which minimize its mean- squared symbol sync jitter. 

A derivation of tbe MAP estimator of symbol sync is given in 
Ref. 1 where an arbitrary time -limited symbol pulse waveshape 
p (t) of energy E = ST and equiprobable symbols are assumed. 
Furthermore, the derivation assumed that the symbol length T 
is known (and remains constant over the obsei'vation interval), 
and that a stored replica of the transmitted signal waveshape is 
available at the receiver. 

To apply this body of theory to the case of Manchester coded 
symbols, we define the pulse shape p (t) by 


pjt) 

s 


VS; O^t^T^/2 

T/2^t^T 
s s 




Also, to allow for average symbol transition densities other than 
one-half, we can assume a data stream of independent symbols 
with a priori probabilities p and q = 1-p and modify the derivation 
accordingly. Without belaboring the details, one can show that, 
analogous to Eq. (9-14) of Ref, 1, the MAP estimate of e-- that is, 

A 

G--is the value of e that maximizes 




K 


A(y, e) =ln f(y, s) = \ In 

k=0 


p exp 


J yk(t)Pg {i- 1 - e]dt 


+ q exp 


'N 


L^o "T, (s) 


J yj^(t)p^[t-(k-l)T^-e]dt 


0 


( 2 ) 


^Many other possibilities exist for statistically modeling the data 
stream; however, the MAP estirrate solution becomes quite 
complex for any model othei* than the one assumed above. 
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wlaere 


K 


y(fc) rs 


aj^PgCfc-(k-l)Ts-el + n(t)j 0 £:kT 


(3) 


k=0 


is the recG ••'•t; ’ symbol data stream in noise over the observation 

interval (0 ^ t ^kT ), is the polarity (±1) of the k transmitted 

^ ^ th 

symbol, and i-S k subinterval defined by (k-l)T^+e^t & 

kT +e , Furthermore, the additive noise, n(t) is assumed to be 
s 

white Gaussian with single -sided noise spectral density 

Alternately, the MAP estimator is that value of e which 
satisfies 

K . (. Sp [t-(k-l)T -e] 


5A(y, e) 
Be 




k=0 




y(t) 


Be 


dt 


+ y{kT_^-fe)p (T ) - y[(k-l)T ielp (0), 


P 


X 


JT^(g)y{fc)^&-(k--l)T^-eldt| -qexpj^ ] 


N 


pexpj“ g.y(t)Pg[Mk-l)T^-G]dtUqex^ . y(t)pJt-(k-l)T^-e]dt | 

t 0 h \ ' [ 0 k' ^ 


= 0 


(4) 


For any estimate of e other than the MAP estimate, the function 
3A(y,e}/Be will be either positive or negative depending on 

A A 

whether e<e or e> s and hence can be used to provide search 
direction, A simple example of a closed-loop symbol synchronizer 
that incorporates BA(y, E)/Be as an error signal is illustrated in 
Pig, 3-1, where for simplicity we have assumed p (0) = p (T) - 0, 
For large F^/Nq, the following approximation can be made* 
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p expj 

\^L (- y(fc)pJt-(fc-l)T^4]dt 
0 k 

;|-q exp- 


P expj 

Ok 

Uq expj 



" sgn]“ J y(t)p„[t-(k-l)T -e]dfc[ 


{5} 


where sgn( ) is the sigmim fiuiction, Likewise, for small E /N^^, 

s u 

we have 


p expj 

li|;L ,~,y(t)pJt-(k-l)T^4]dt| 
O k ^ 

-qexpj 


1 

pexp- 

1 


:|+q exp 



- (p-q) + 4pq 


^ f ^ y(t)p [t-{k-l)T -e]dt 

L^O ® " 


(6) 


Using (B) and (6) in Fig. 3-1, gives the equivalent high and low 

E^/Nq closed-loop configurations illustrated in Fig. 3-2. 

For the Manchester pulse of (1), the output of the minus one 

amplifier following the differentiator in Figs, 3-1 or 3-2 is illustrated 

in Fig, 3-3, The strength I] of the impulses at the edges of the 

symbol pulse is, in general, a random variable and depends upon 

the polarity of the preceding and succeeding symbol pulses. For 
til 

example, in the k symbol interval, Tj = 1 if 1 ” ^k ~ ^k+1^ 
and 0 if a^^ ^ ^ and ^ Since the probability of a 

data symbol transition is p = 2pq and the probability of no 

li 

transition is 1-p,, then the average value of T], namely T, is 
given by 


~ 0 XP<. + 1 X (1“P.) = ^ - Pf 


(7) 
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Fig, 3-3, Manchester Syrnbol Pulse Shape and Its Derivative 


— /I .. . 

Note that Eq, (7) is valid independent of the value of E /N„. 

In practice, the impulse ftinctions discussed above would be 

approximated by finite width pulses. With regard to the impulses 

at the edges of the symbol pulse, there are two natural ways of 

making this approximation each of which leads to a different 

closed loop symbol synchronizer configuration. In one case, the 

cross -correlation reference signal at the lower phase detector in 

Fig. 3-1 or 2 might appear as the waveform illustrated in Fig, 

3 -4a. This reference signal has the effect of gating the integrate- 

and-dump on for /2 sec at the beginning and end of the symbol 

s 

interval and twice as long at the mid -symbol position. According 
to MAP estimation theory, the fractional "window width, " 
should he chosen as small as possible^ Frora a practical stand- 
point, however, we shall soon see that other considerations are 
involved in the selection of its value. In the second case^ the 
same, cross -correlation reference signal would appear as in Fig. 

3-4b. Here, the edge gates overlap two adjacent symbols. 

Regardless of which impulse function approximation is adopted, 
in a practical symbol synchronizer implementation, the I'unning 
accumulator could be replaced by a digital filter whose memory 
is chosen relative to the value of K (Ref. 6), Whereas, the input 
to the humped phase oscillator was a linear combination of the 
previous K+1 outputs of the imiltiplier, we now have an accumula- 
tion of the infinite past of these quantities each weighted hy the 
impulse response of the digital filter. By adjusting the memory 
of this filter, the input to the humped phase oscillator can be made 
to essentially reflect only the previous K-f-1 multiplier outputs.. 

Since in practice the signal epoch will never remain Constant in 
time, a weighted sum of the data that accentuates the recent 
symbols and attenuates the others should be more appropriate than a 


■REPTf.r>T)TTnmTT,TW Q T? 

ORIGINAL PAGE IS POOR 


- 47 * 








On the basis of the foregoing information, we propose to study the 
steady-state tracking performance of the symbol synchronizer 
implementations in Fig. 3-5, Those in Figs, 3 -5a and 3-5b are based 
on the approsrimation of Fig. 3-4a while those in Figs. 3-5c and 3-5d 
make use of the approximation of Fig, 3-4b, The tecliniques 
which will be employed in carrying out this study are analogous 
to those given in Ref- 1 where they are used to characterize the 
performance of symbol synchronisers with NRZ inputs. Namely, 
under the assumptions that the input timing offset G is essentially 
constant over a large number of symbols and that the loop response 
is very slow with respect to a r'i,‘mbol interval (W T *=5^ 1 where 
W is the two-sided loop bandwidth), the symbol synchronizers of 
Fig, 3-4 can be modeled as the continuous phase-locked loop 
given in Fig. 3-6. Developing the equivalence between Fig. 3-5 
and Fig. 3-6 relies on finding 1) the average loop S-curve g(X ) as 

A 

a function of the normalized timing X = (e-e)/T and 2) the two- 

s 

sided spectral density S((ju, X) of the equivalent additive noise n (t), 

A. 

Once having determined these quantities, the steady-state 
performance can be found by application of the Fokker-Planck 
equation. 

The approach taken in finding g(X) and S(UJ, X) follows that 
given in Refs, 1, 7, and 10, In particular, if one replaces the 
samplers in Fig, 3-5 by sample -and -hold circuits, then the 
input to the digital filter becomes the continuous stair case-t}'pe 
waveform e(t) defined by . 
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^ A 

e(t) = kT+e^t ^(k+l)T + e 


( 8 ) 


whose statistical mean cohditioned on fixed X is g(X) and whose 
power spectral density conditioned on fixed X is S{tu, X). Mhthe- 
matically speaking. 


gM 

n^{t) 


= e(t)-g{xr 


S{tu,X) = ^<R{t,T;X)>} 


R(t, T;X) ^ Jnj^(t)nj^(tf T) lx} (9) 


■ w^^here E^^ ^ represents the conditional expectation on X both with 
respect to the noise and the signal (symbol sequence), the symbol 
3’ denotes the Fourier transform operation, and the symbol y 
denotes time average. Furthermore, to evaluate the autocorrela- 
tion <CR{t, TjX)>, it is sufficient to determine the discrete function 
RQ(m, X); m = 0, ±1, 12, , , , defined by 

(10 

Then, since n^(t) is conditionally cyclostationary (i. e. its auto- 
correlation function R(t, t;X) is, for fixed T and X, periodic in t 
with period T^, we have 

R.(m+l,X)-R (m, X) 

<^R(t, t;X)}> = Rp(m, X) +— — — — ^ — (T-mX^); (11 

s 

mT ^T^(m+1)T 

... s s . 

m = 0, ±1, ±2, ... 

The autocorrelation function RQ(m, X) has the following properties: 
(1) R^(m, X) has nonzero value only at m = 0, ±1, This is a 
consequence of the independence assumption made 
regarding the data symbols. 




( 2 ) 

Thus, Eq, (9) combined with (10) simplifies to 

2T 

S(iu, X) “ ^ <CR(^:i T;X)> cos U)TdT 


r 


R|^(0,X) + T COS OJTdT 

s 


2T r R^^(1,X) -] 

■1 2j |R|^{1, X) - (t-T^) cos 0)TdT (12 

*^s ^ ® ■ 


if ® ] 

L 


Since, as before, it is assumed that T 1, it is sufficient 

3Li s 

then to consider only the value of S(U), X) at aero frequency --that 
is, S{0, X)--and assume a flat spectrum of this value for all U) of 
interest. From (12), we get 

S(0,X) tJRq(0,X) + 2R^(1,X)1 (13 

Furthex'more, since in any practical system little data degradation 
due to imperfect symbol sync can be tolerated, then extreme 
accuracy is required in establishing symbol synchronization. 

Since this implies very large loop signal-to-noise ratios, then 
the value of S(0, X) is essentially the noise spectral density seen 
by the loop at X == 0, viK, , S(0, 0), see Ref* 1, Gh. 9. 

3.4 Application of the Fokker -Planck Technique to Obtain 


Performance 

The stochastic differential equation which describes the loop 
illustrated in Fig, 3-6 is 

i = -K^F(p)[g(>.) + (14 

where p is the Heaxuside operation and the dot denotes differentia' 
ti on with respect to time* In writing the above, we have assumed 











zero static timing error which is consistent with our previous 
assumption that e(t) 0, 

In dealing with equations of the above form one can, at best, 

hope to find the probability distribution of X(t], viz,, p(X, t) from 

2 

which the mean-squared timing jitter cr^(t) can be calculated. 

Under the assumptions previously made regarding the equivalent 
additive noise process n.(t), the solution for p(X, t) can be found 

A 

via the Fokker-Planck method. In particular, for the first-order 
loop case, i. e, F(p) = 1, X(t) is a first-order Markov process 
whose probability density function satisfies the diffusion equationj 
Ref. 1, " 

Sp(X.,t) , ,2 

~St ajT ^ ^ (15) 

0 

where t) and K^(Xq, t) are defined by the conditional expecta- 

tions 

K (X ,t) = lim ^ E[[X(t+At)-X{t)]|X(t)=Xj (16) 

^ At^O 

K-(X^,t) )= lim ^ E{[X(t+At)-X(t)]^lX(t)=X } (17) 

Of interest is the case where p(X, t) converges to with time to a 
stationary probability density function independent of the initial 
condition, which can be used to evaluate the steady- state behavior 
of the symbol synchronizer. Denoting 

p(X) = lim p(X, t) (18) 

■ ^ca 

CTn terms of the time discrete symbol synchronizer configurations 
of Fig. 3-5, the first-order loop assumption is equivalent to assum- 
ing a digital filter with unity weighting coefficients. 







we obtaiin froirt (15) tbe stationary equation 


^[Kj{X)p(X)] = [k^WpW] 


where from (9), (13), and (14), we obtain 


K^(M = -K.^g(M 

K-{\) = K^S(0,X ) = K^S(0, 0) 

“V V 


Substituting (20) into (19), and solving, gives the probability 
density function 


p(X) = exp 


K^S{0, 0) 


; 0 


P(-X) = P(X) 


where Cj^ is a normalization constant chosen such that 


p(X)dX = 1 


Now for large group signal- to -noise ratios, g(X) can be 
approximated by 


g(X) = 


^/SK X 
g 

STK X 
g 

^SKX 


STK X 
g 


(Fig, 3-5a) 
(Fig. 3 -5b) 
(Fig. 3-5c) 
(Fig. 3-5d) 


where K is the slox>e of the normalized S-curve, g (X), at the 
§ . ^ 

origin. Then defining the ratio of symbol rate = 1/T to 

S . S'"' 
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single- sided loop bandwidth. 


W. 


we have for large R. 6 (bnt necessarily large R ^ that 

^ cs o ■ S 


^ f \%(X)dX^ 
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(Fig, 3 -5a) 


(Fig, 3-5b) 


(Fig. 3 -5c) 


(Fig, 3-5d) 


The value of symbol edge gate amplitude y (see Fig. 3-5) may 

now be selected to minimize ct^. Differentiating (25) with 

respect to Y an<i equating to zero gives 

A-2pq (Fig. 3-6a) 

l-2pq (Fig. 3-5b) 

2{l-2pq)-|(2pq)V(R )7TT‘e3cp(-R J/erf^S:? 

S - o ' 5 


U5R^(2pq)[l-erf -v^l-2?(2pq)^“^ exp(-R^)erfv^--^ exp (-R^) 

(Fig. 3-5c) 

(26c 


2(l-t2R )[(l-2pq) - i5(2pq)l 

' S • 

l+4(l-2pq)R -:^?[l+4{2pq)R ” 
s . s 


(Fig. 3-5d) 
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Kofce that the value of Y which mlmtnizes for Fig, 3 -5a and 

Fip;, 3 -5b is independent of R and is identical to the value of T] 

s 

[seeEq. (7)] suggested by the MAP estimation theory. The 

values of Y . lor Figs, 3-5c and 3-5d are, in general, functions 
opt 

of R , ? and 2pq, Fig, 3-7 illustrates y of Eq. (26c) vs R for 
s ^ opt s 

fixed I and p - 2pq as a parameter. Keeping in mind that Fig. 

3-5c is intended for use at large signal- to -noise ratios, we note 

from Fig, 3-7 that over a wide range of R , Y . can well be 

s opt 

approximated by its value at R = i, e, 

s 


opt 


= Z(l-2pq) 


|R s=co 
S 


Again, note the similarity with the result predicted by the MAP 

estimation theory. Figo 3-8 is a similar plot to Fig, 3-7 for 

Y pj. given by {26d). Here we see that over a similar raiige 

of small signal-to-noise ratios, is well approximated by its 

value at R = 0, i, e» 
s 


opt 


R =0 
s 


2r(l-2pq)-lg(2pq)l 
" " 1-1 


(27) 


(28) 


From (25) and (26) we see that a comparison of the various 

symbf*'! synchronizer configurations on the basis of minimum 

2 

mean- squared timing Jitter can be achieved by studying the 

behavior of the two ratios 


2 

mm 

'“x . ’c 

min 


2pq) 

1 + (l-2pq)^ 


(29) 


and 


<“x . >b 

min 

»; . I. 

mm 


{5.R„>2pq) 

. J-i ° 

1 + (l-2pq)^ 


( 30 ) 
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where 


- 1 + 


(31a) 


. [(l-2pq)erf^-i5(2pq)j ““ exp(-R^)']' 


TK 


er/v^{i+|?Rg(2pq)[l-erf^V®^]+?(Epq)^-:~ exp(-R^)er£yB? exp (-R^)} 

,2 


^ 2(U2R^)[{i-2pq)-i5(2pq)7 

^ 1 -f 4{l-2pq)R^ + ^5 [1 + ^(2pq)R~T 


(31b) 


Figs. 3"9 and 3-10 illustrate these two ratios (in dB) as a fiinction 
of R with § fixed and p = 2pq as a parameter. "We observe from 

S t 

Fig. 3-9 that over a large range of values for R^, and p^, the 
minimum mean-squared timing jitter associated with tie synchro- 
nizer implementation of Fig, 3-5a is larger than that associated with 
Fig. 3 -5c. In fact, over a wide range of large signal- to-noise ratios, 
the ratio given in (29) can be approximated by' 


K L 2 

min 1 + 2(l-2 pq) 

2 


(32) 


K : >c 

min 


1 -t- (l-2pq) 


Similarly, from Fig. 3-10, we observe that over the range of 
parameter values considered, the symbol synchronizer configura- 
tion of Fig* 3-5d outperforms that of Fig, 3- 5b, 

To complete the picture, we compare the mean- squared jitter 
performance of Fig, 3-5a with Fig. 3-5b and Fig, 3-5c with Fig. 
3-5d so as to determine their relative behavior in the region 
between small and large signal- to -noise ratios. The motivation 
behi,nd such a comparison is the fact that in this region of medium 
signal-to -noise ratios, the appropriate nonlinearity to be used at 
the output of the in-phase arm integrate-and-dump circuit as 
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Fig. 3-11. A Compai'isoii of the Mean-Squared Timing Jitter Performance of Figs. 3-Sa and 3*--.h 







suggested by tbe MAP estima.tion theory is neither the sigrmm. 
nor the linear function [see Eq. (4)] . Thus, it is not immediately 
evident which implementation is preferable in this in-between 

signal- to-noise ratio region. Pigs, 3-11 and 3-12 illustrate,. 

Z 2 2 ,2 

respectively, the ratios {a. ) and {a K 

A.o A., a A. a. A.C 

„ e -I min^ mm min mm 

vs, R With 5 and p as parameters, 
s c 

3.5 Further Numerical Results 

The S-curves (phase detector characteristics] for the four 
symbol synchronizer implementations of Fig. 3-5 are illustrated 
in Figs. 3-13 to 3-1 6 respectively for p = . 5, various values of 
window width §, and R^ as a parameter. We observe from these 
illustrations that over a certain range of signal-to-noise ratios, 
the various loops exhibit a stable lock point which is one half 
symbol (X= .5) away from the desired lock point at X = 0. 

Although, this midsymbol lock point is less probable than the 
one at X" 0, an attempt should be made to resolve it whenever 
possible. For sufficiently large signal-to-noise ratios, viz.j 
greater than 10 dB, the midsymbol lock point vanishes for the 
synchronizer implementations of Figs, 3-Sa and 3 -5c. On the 
other hand, the S-curves for Figs. 3 -5b and 3-5d are (for 
p. == . 5) independent of R . For smaller and smaller transition 

L S ' 

densities, p^, the S-curves become more and more symmetrical and 
in the limit of zero transition density, i,e. a square-w;ave data 
waveform, the S-curves are periodic with period equal to 1/2, 
and thus, the midsymbol and desired lock points are equiprohable. 
Therefore, some form of ambiguity resolution circuit Will be 
required with the use of these phase detectors, 

3.6 Cycle Slipping Performance 

Up until now, we have treated the normalized timing error 
process, X(t), in the symbol synchroniser as a stationary 

Ph 

a^Lin 


i^om. 
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Fig, 3-l3a, S-G-tirves for Symbol Synchronizer of Fig, 3-5a; p^. = . 5, S = , 5, Y = Yopt 
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process possessing a bounded variance [see Eq, (25)] , This 

resulted from the fact that we have in reality been discussing 

the behavior of the timing error reduced modulo T , or 

s 

equivalently, the normalized timing error reduced modulo 1, " 

completely describe the \(t) process, one must account for the 
component of its variance that results from diffusion, that is, 
cycle slipping. The steady -state effect of cycle slipping is 
perhaps best described by evaluating the diffusion coefficient, 
that is, the rate at which the variance of the timing error is 
approaching infinity, or equivalently, the total average rate of 
cycle slips, S. The cycle slipping phenomenon in symbol 
synchronizers is not xinlike the same phenomenon in phase - 
locked loops. In fact, with appropriate normalizations, the 
expressions for S given in Ref, 1 for the phase-locked loop can 
be easily generated to apply to the symbol synchronizers of 
Fig. 3-5, In particular, it is straightforward to show that 


where 


A 

V U ^ 

^nO^^^ “ K 


n ^ — 

P " 2 


with Q determined from (25), The expressions for the loop 

/ L__ . . • ■ , 

S-curves given in (14) can be integrated in closed form; however 
the remaining integrals in (33) must be performed numerically 
on a digital computer, ‘When this is done, one finds that for all 
values of R S . 1, 5 a 100, §^.5, and a . 05, the normalized 






average rate of cycle slipping as computed from (33) is below 
lO Tbus, for most practical applications, cycle slipping 

will not be a problem, 

3, 7 Symbol Synchronizer Performance in the Presence of Clock 



In deriving the results given thus far, we have made the 
assumption that the unknown epoch, e{t), of the input symbol 
data stream is constant with time. Quite often in practice, e(t) 
0 Qi-eit/T. Under these conditions, the stochastic differential 
equation of (IB) becomes 

j. = - K^F(p)[g(J.) + I 

and the symbol synchronizer operates in the presence of a clock 
frequency offset. Applying the Fokker -Planck technique to (36 ), 
we obtain a generalization of ( 21 ), namely, 

r 2 jg(X)dX" 

p{h) = exp 0)' “ K^S(0, 0) ' 

where G_ is a normalization constant. For large R 6 , the 
mean-squared timing jitter is still given by (25); however, the 
stress in the loop caused by has, as we shall see shortly, 
a pronounced effect on symbol error probability performance. 
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4.0 AVERAGE ERROR PROBABILITY PERFORMANCE OF 

convolutionally coded systems in the presence 

OF IMPERFECT SYMBOL SYNC 


4.1 Introduction 


Digital communication systems whick, by necessity, require 
symbol syncbronization for successful operation, suffer a perform- 
ance degradation due to tbe timing jitter associated witb tbe symbol 
syncbronizer. Accurate assessment of this degradation is essential 
in predicting overall end-to-end system bit error rate performance. 
Regardless of whether the system is coded or uncoded, proper 
design of the symbol synchronizer loop plays an important role in 
reducing the degradation to a minimum. For a given symbol 
signal-to-noise ratio, R , symbol rate /? = 1/T , and loop band- 


width W , it is desirable to select a configuration which yields 
L 


minimum mean- squared timing jitter and is commensurate with 
the requirements on irnplementation complexity. Several such 
configurations for Manchester coded symbols were described in 
the previous section. 

"When the system is uncoded, it is sufficient to know the 
probability density function, p(X}, of the normalized timing error 
\ in order to analytically evaluate the error probability perform- 
ance degradation due to a noisy symbol sync reference. In 
particular, assuming that the timing error is essentially constant 


over the symbol duration T^^, then the bit error probability P^ 


is given by 


P {b)p( ;dX 
1 r, , 

3 


where 


Pj^(X) - Ir, e: fct,/ir(1.2ixrj -i -‘(.1 .pj^)errc[v'pr(l-4|x|] 


refroducbii tb ■' ^ 
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I 


and p is again the symbol transition density, 

U 

For convolutionally coded systems, the effect of imperfect 
symbol sync on bit error probability performance is, in general, 
more difficult to evaluate analytically and, indeed, in many cases 
must be determined by either simulation methods or hardware 
measurements. Fortunately, however, for the Shuttle S-hand 
relay links, suitable assumptions regarding the behavior of the 
system can he invoked and an approximate analytical evaluation is 
possible. In the next part of this section, we describe such a 
technique and present numerical results for the Viterbi decoded 
K=7 rate 1 /3 Shuttle code, 

4.2 An Analytical Evaluation of Performance Degradation due to a 

Noisy Symbol Sync Reference 

Generally speaking, the presence of a timing error in the 
symbol synchronizer causes a reduction ill the signal-to-noise 
ratio available per symbol at the Viterbi decoder input. If the 
loop bandwidth is large relative to the symbol rate, so that the 
timing error process varies rapidly over the duration of a single 
symbol, then the degradation in average symbol signal-to-noise 
ratio is a good measure for estimating bit error rate degradation 
in Viterbi decoding. Unfortunately, in rro st practical convolutionally 
coded systems, the loop bandwidth of the symbol syiichroniaer is 
considerably narrower (perhaps hy several order of magnitude) 
than the symbol rate, and thus a more direct approach, namely 
direct computation of the bit error rate, is required. 

If we examine the symbol synchronizer configurations 
previously discussed, we observe that the output of the in-phase 
arm integrater-and-dump when passed through a soft quantizer 
produces the soft decisions required at the Viterbi decoder input. 


p. 




For siinplicity of presentation, we shall assume an infinite bit 
quantizer and, thus, the integrate and dump output is fed directly 
into the Viterbi decoder. Using previously defined terminology 

A 

the integrate -and -dump output z(t) at time IcT is given by 

s 

z = z(kX +s) = Tp I ^ y(fc)p 0:-(l^-l)T-e]dt 

s^k-l)Tg+e 

( ^[a,(l-3|xl)-a_, |x|] + n.;X >0 




where 


^ kT +e ^ 

ru i — f ® nWp [t-(k-l)T-e]dt 

■'(k-l)T +e 


XN 

E{nj^l = 0; i E{n^) - ~ 


getting k =ivtm (1/v is the code rate), then a. specifically 
th 

represents the m (m = 1, 2, . , *m) encoder output symbol 

fch 

cor X'esponding to transmission of the ^ information bit. 

The Viterbi decoder, being a maximum -likelihood decoder, 
computes the test statistic (metric) 


L r V 


'I " > . ^k^k 


z a 
Xv+m 


JL^l ^m=l 


for all possible sequences of encoder output symbols corresponding 
to the transmission of L information bits and chooses the input 
message that yields the largest value of q. The brackets in (6) 






indicate that the computation of q '.a done or? ».ree branch at a 
time as one sequentially steps through the L branches of the tree. 
By choosing the value of B to be several times the constraint 
length of the code, the selected path through the tree can be used 
to decode the first bit in the sequence of L, 

Using the decoding metric of Eq. (6), a tight upper (union) 
bound to the average bit error probability is (Refs, 1,2) 




^ > n.m.P. 

j-d , 
mm 


( 7 ) 


where P, denotes the probability of error in comparing the 
3 

correct path to an incorrect one which differs in j symbols from 

it, d . is the minimum free distance of the code, m. is the munher 
min j 

of paths at distance j (differ in j symbols from the correct path) 

and n, is the number of bits each of these paths differs from the 
3 

bit sequence corresponding to the cori’ect path. In the absence of 

symbol sync error, P. is given by 

3 


P, = erfc 


( 8 ) 


for each of the m. paths at distance j. When symbol sync timing 
3 

error is present, the specific location of the j symbols which are 
in error along the Incorrect path at distance j also affects the 

calculation of P.. Thus, P. is, in general, different for each of 

3 . 3 . . 

the m. different incorrect paths at distance j.‘ To see this, let 
fil ^ 

{aP ] , k = 1, 2, . . , , L, i = 1, 2, . . . , m. denote the sequence of 

th ^ , 

symbols along the i incorrect path at distance j. Such an 

incorrect path will be chosen if and only if its metric is larger 

than that of the correct path, i. e. 






or 


Lv ... IjV 


z 3. 
k Ic 


( 9 ) 


k=l 

IvV 


k=l 


k=l 




( 10 ) 


fi) 

Let S, (j) further denote the set of locations where ^ a^* 

Clearly S.{j) contains j elements, Tlms> from (10), we get that 
th ^ 

the i incorrect path will be chosen if 


or equivalently. 


„E 


z, a. >0 
k k 


(IM 


keS.(j) 


z^ a, <0 
k k 


( 12 ) 


keS Jj) 


Substituting (3) into (12), we find that the probability of selecting 
th 

the i incorrect path at distance j from the correct path is 


P.(i) = Pr 
3 


keSJj) 


Vkdi 


Wl< 


n. a 
k k 


kGS.(j) 


(13) 


n a, is Gaussian with 
k k 


where the ±is determined by the sign of X in accordance with (3). 
Since, from (5) the random variable 
zero mean and variance jNp/2T, we i get that 

— 


P (i) = I erfc 

J 


Finally, then, for a given value of normalized timing error 
the upper bound on bit error probability analogous to Eq. (7) is 


(14) 


I pti PimPPCIBILfg 

ORIGINAL 
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p^(X) ^ 






(15) 


j=d . i=l 
min 


The summation in (14) is upper and lower bounded by 


keS.(j) 

and for reasonably large j is well approximated by 


(16) 


keS.(j) 


■ [l“3lx|-a^aj^^|xl] = j[l-3lxl-W(i-PfeHUIPfc] 

== j[l-4lx|+2]X|pj.] 


(17) 


with the approximation improving as j increases. Since in comput- 


ing the upper bound in (15), the minimum value of j 


since d . =14 for the K=7, \^=3 Shuttle code (Ref, 3), we shall 

mm 


use (1?) in (14). When this is done P.(i) becomes independent of i 

3 


and (15) simplifies to 


C3 


P^(X) 5 


n.m (I erfc{ [l-4|x| 4 2 |xl p ]}) (18) 

J J S 1. 


•j=d . 
min 


which is identical to (7) (after substitution of (8)) except that JK ' 

■ s 

is replaced by (l-4-lX|+2 jxjp.)* The average bit error 
s u ■ 

probability is then uppei' bounded by 


09 


P^ ^ J n.m^ I erfc{,v^^[l-4|x|+2 j x[pj.]]p(X)dX (19) 

^ Ji, J J 


3=d . 

mm 


Since the union bound of (7) is very tight over a relatively wide 
range of bit signal-to -noise ratios R^ = vR^ at low decoding error 
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rates, the polynomial curve fit to fclie simulation results may be 
used in place of of (18) by making in it the substitution 

//R^ {i-4|?tj + 1 mp^) for VR|^. Finally, then, 

1 

S f ^ P{,(X)p(Ma>- (ZO) 


where 


log^^jP^(X) = -.1206.(10 log^pRyX)) ~.50(l01og^QR^(\)) 


-1.0706 


and 


R^(X) i R^(l-4lXl + 2lxlp^}^ 


( 21 ) 


(22) 


Using the symbol synchroniaer configuration in Fig, 3-5a, 

Figs, 4-1 to 4-5 illustrate P^, as determined by the numerical 
integration of (20), vs R for p = , 5, 5= .25, and values of 

D U 

normaliaed clock frequency offset 6f/B = 0, , 02, , 04, . 06, and , OS. 
The ratio of symbol rate to loop ba’^dwidth is a parameter on each 
of these curves. The limiting curves for 6 -4» are determined 
from (21) with the substitution X = Af/4B^. For all practical 
purposes, this limiting performance is obtained with a 6 of 2000, 
Figs. 4-6 and 4-7 illustrate the bit error probability performance 
for other values of transition density and no clock frequency offset. 
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Fig. 4-7. 
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Bit Error Probability Perfonnance of lC-7, Rate 1/3 Shuttle Eode 
in the presence of Noisy Symbol Sync; ^= .^5 u = .S^ 'kf/B » 0 
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